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SATELLITE-TRACKING AND EARTH-DYNAMICS 
RESEARCH PROGRAMS 
Semiannual Progress Report No, 29 

1. INTRODUCTION 

Tills report covers the following activities in Smithsonian Astrophysical Observa- 
tory's (SAO) earth-dynamics programs; 

Ac Satellite- tracking network operations (NASA/OTDA) (see Section 2). 

B. . Satellite geodesy and geophysics programs (NASA/OA) (Section 3). 

C. Atmospheric research (NASA/OSS) (Section 4). 

The following paragraphs outline some of the highlights of our program activities 
during the reporting period. 

During the last 6 months of 1973, approximately 46, 000 successful range measure- 
ments were acquired by the SAO laser stations in Peru, South Africa, Brazil, and 
Arizona; this constitutes a 50% increase over the comparable period of 1972, 

The Peole satellite-tracldng campaign conducted in conjunction with the Centre 
National d'Etudes Spatiales (CNES) was completed in August 1973. Tracking of Peole 
began in April and continued after the D5B launch vehicle failed. During the campaign, 
the SAO network obtained 4482 validated returns of 310 arcs of Peole, This is an 
impressive figure, compared to only 1178 Peole returns on 195 arcs acquired during 
the International SatelUte Geodesy Experiment (ISAGEX) and the Earth Physics 
Satellite Observation <’ampaign (EPSOC), This improvement is attributed to new 
prediction software and operating techniques. These data are of particular value 
for obtaining more accurate gravity-field and zonal-harmonics coefficients. 
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Work continued on tlie design and early construction phase of the prototype pulse 
digitizing system, which presents, in digital form, the outgoing and return pulses 
from each laser tracking system. We plan to install a pulse digitizing system at each 
SAO laser station in order to improve ranging accuracy in preparation for Geos C. 

The Analytical Satellite Geophysicc Department began work on the Research and 
Technology Operating Plans (RTOPs) assigned to SAO in support of the Earth and 
Ocean Physics Applications program (EOPAP). These are discussed more fully in 
Section 3. 

The first half of Fiscal Year 1974 constituted a transition from efforts supporting 
the National Geodetic Satellite Program (NGSP) to those for EOPAP. In addition, a 
number of NGSP activities wore completed. SAO’s contribution to the NGSP, which 
also comprises the 1973 Smithsonian Standard Earth {HD (SE III) (Gaposchkin, 1973), 
was submitted to the American Geophysical Union and will be published as Chapter 9 in 
the NGSP Pinal Report. A comparative study of all NGSP contributions is now under 
way. The preliminary results of this comparison indicate that geocentric coordinates 
have quoted accuracies that are optimistic by factors of 2 to 5. Even so. Independent 
estimates of coordinates agree with each other to within 10 m, the NGSP goal. 

The results obtained from work on some of the RTOPs are being used as support- 
ing material in preparing SAO’s proposal to OTDA for upgrading SAO's laser systems. 

SAO participated in a number of Lageos activities during this period. This 
included a review meeting on Lageos program objectives and scientific and technical 
features at National Aeronautics and Space Administration (NASA) Headquarters on 
11 October, followed by subsequent meetings both at NASA Headquarters and at 
Marshall Space Flight Center (MSFC) in Huntsville, Alabama, Discussions covered 
various views concerning the orbit altitude, orbit inclination, and the satellite size and 
mass required for the specification of the Lageos system. From these data, NASA will 
compile the Phase B definition study. A number of memoranda were provided to 
NASA/OA as a result of analyses performed by SAO - e. g. , one by Dr. E. M. 
Gaposchkin entitled ’’LAGEOS Orbit Specifications, ” dated 10 December 1973, and 
others on the Lageos retroreflector-array requirements, by Dr. M. R. Pearlman 
and Mr. D. Arnold. 
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Dr. L. G. Jacchia's attempts to model the variations in thermospheric composition 
as observed by satellite mass spectrometers, and to reconcile these variations with 
doppler temperatui’e data and total densities from satellite drag, have scored remark- 
able success. T^vo simple formulas have been found to represent the observed seasonal- 
latitudinal variations in composition and the diuimal variations of the individual atmos- 
pheric constituents at any given location and lieight. Each equation contains a single 
numerical constant. Attempts are also being made to find a way to represent the global 
variations of temperature and composition with geomagnetic activity. 

A study of the effect of radiation pressure and atmosplieric drag on the balloon 
satellite 19G3 SOD has been completed. Mr. J. Slowey has found that by using a single 
value of the reflection parameter, it is possible to model perturbations due to direct 
radiation pressure in all the elements. 

Slowey has be^ an work on the determination of atmospheric rotation from the 
secular decrease in the inclination of satellite orbits. Most of the necessary program- 
ing has been completed, and the first results were being obtained at the end of the 
reporting period. 
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SATELLITE-TRACKING NETWORK OPERATIONS (STADAD) 


2, 1 Satellite Obsorvlng Campaipts 

SAO contimies to track Geos 1, Geos 2, and BE-C in support of its ongoing NASA- 
supported research programs in geophysics. The data acquired during the last 2 years 
on these satellites are now being analyzed by SAO's Analytical Satellite Geophysics 
Department to measure the motion of the earth’s pole (Chandler motion) and to refine 
geodetic parameters in preparation for Geos C. 

The Peolc satellite-h’acking campaign, initiated 11 April 1973 in conjunction with 
the CNES, was completed on 7 August 1973. Tracking of Peole began during the pre- 
launch phase of a proposed D5B tracking program and continued after the failure of 
the D5B launch vehicle. During the campaign, 4482 validated returns on 310 arcs of 
Peole were obtained by the SAG network. In comparison, only 1178 Peole returns on 
195 arcs were acquired during ISAGEX and EPSOC. This improvement is attributed 
to new prediction software and new operating teclmiques. These data are of particular 
interest for obtaining more accurate gravity -field and zonal-harmonics coefficients. 

2. 2 Laser Data 

During the reporting period, nearly 46, 000 successful range measurements were 
acquired by the SAO laser stations in Peru, South Africa, Brazil, and Mt. Hopkins 
(see Table 1). This represents a 50% increase over the last 6 months of 1972. The 
Geos 1, Geos 2, and BE-C satellites were tracked 7 days a week until September, 
when the laser schedule was reduced to 6 days for operational reasons. 

Afteiv their summer rainy season, the Toliyo Astronomical Observatory (TAO) began 
operating its laser system again in November with a modified prediction. format. An 
increase in data has been seen. 




BlANIv biux 
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Table 1. Successful laser measurements (numbers in parentheses indicate successful arcs). 
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A now coude lasex’ mount wns puvcliasod by the National Technical University (NTU) 
in Athens. It wns temporarily sot up for testing at the former BaUor-Nunn site wlille 
the laser building was being modified to accommodate the coude mount. In October, the 
laser and electronics wore sent to Othoni, an island northwest of Corfu, to obtain range 
measurements between Corfu, Othoni, and Italy. The purpose of those moasux*emcnts 
was to obtain a more accurate geodetic tic between Greece and Italy and, In turn, to 
estublisli a more accurate geodetic location for Gi’oeco in the Buropean datum. Pre- 
liminary results indicate that the project was a success. 

The laser system in Ethiopia, on loan from ONES, was not operational during the 
reporting period. Persistent electronics and supply problems liave plagued that sys- 
tem. At a meeting held with a representative of CNES in December, it was decided 
that the system would not be operated without a commitnxent of engineering support 
from CNES. SAG has now suspended the laser operations activity in Ethiopia with the 
understanding that the French may provide support for special observing periods in 
Calendar Year 1374. 

2. 3 Baker-Nunn Camera 


Ten Baker-Nunn cameras operated by SAG and its cooperating agencies made over 
8000 successful observations of 22 satellites (see Table 2). The cameras provided 
routine observations to support the generation of pointing predictions for the lasers 
belonging to SAG and its cooperating agencies. 

Table 2. Total number of Baker-Nunn observations by month. 


Month 

Number of 
observations 

July 

1561 

August 

1542 

September 

1539 

Getober 

1445 

November 

982 

December 

934 

Total 

8003 
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Tlie network made over 2400 successful observations of the five atmospheric 
satellites in support of NASA-sponsored research by Jacchia (see Table 3). Special 
support was also provid<^?H by tracking the Cosmos 382 rocket body for geodetic 
research for Mr. Carl A. Wagner of Goddard Space Flight Center (GSFC). We 
stopped tracldng the object on 15 December at the request of NASA, Tracking support 
condnued on 1971 54A for Dr. Desmond Klng-Helo of tlio Royal Aircraft Establishment 
for his study of 15th-ordor resoriance. In addition, special photographic support on 
satellite 1972 37A was provided for the Air Force. 

Table 3. Baker-Nunn data on atmospheric satellites, 

July to Decern iDer 1973. 


Satellite 

Number of 
observations 

Vanguard 2 

350 

Explorer 8 

445 

Explorer 19 

703 

Ex-plorer 32 

323 

Explorer 39 

662 

Total 

2483 


The Baker-Nunn camera in Australia remains in storage, awaiting construction 
of a building at a new site in Orroral Valley, 

The twelfth Baker-Nunn camera, operated by the CNES in Upper Volta, was dis- 
mantled in October, and the mirror was sent to be resilvered in preparation for the 
launch of the French satellite Starlette. 

The network began nightly observations of Comet Kohoutek on 20 November, and 
by the end of December, 78 successful photographic sequences had been obtained by 
eight stations (see Table 4). 
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Table 4, Scqucncos ol! Comet Kohoutok made from 
21 November to 14 December 1973. 



Si)ain 

India 

Hawaii 

Mt. H. pkins 

Peru 

Novembo r 

4 

0 

0 

3 

2 

December 

_7 

A 

3 

JJL 

G 


11 

10 

3 

14 

8 


South Africa 

Ethiopia 

Brazil 

Total 


Novemljer 

2 

3 

0 

20 


December 

G 

bB 

G 




8 

18 

G 

78 



2.4 Englneeriag 

The second half of 1973 has been devoted to the design and early construction phase 
of a prototype system to digitize the outgoing and return pulses of the laser tracking 
system. This system is being built to improve the laser ranging accuracy in prepara- 
tion for Geos C, 

Various system components, such as the digitizer, paper-tape perforator, and 
pulse processing devices, were selected and purchased during the summer. In the 
fall, SAG built an intercoupler to handle the system flow and to format the data from 
the digitizer, the laser system clock, and the range counter. 

The intercoupler was operational by December, and system integration has already 
begun. Tlie entire system will undergo shakedown and calibration tests early in 1974, 
after wliich it will be installed at Mt. Hopkins for field testing. 

Work continues on the software development for the minicomputers for the field 
stations. The full prediction progx*am that runs from orbital elements is now nearly 
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complelc. Comparisons with pi’ccUcttons generated on the CDC 0400 computer show no 
differences, Present efforts are being applied to the Input/output software and to other 
npe rational support routines (see Section 3). 

2 . 5 Comnuin icatlons 


Experiments with direct communications between SAG and Mt. Hopkins wore 
conducted In November by using acoustic couplers with the FTS phone Ihies. Initial 
problems were solved, and the direct link has boon established, thus relieving the 
Mt. Ilopldns staff of a 20-mIlo drive to retrieve laser and Bakor-Nunn predictions. 

A dh’Cet link was established between STADAD offices and the SAG Commimlcatlons 
Center to speed passage of administrative traffic. 

The rodiv* TTY link with Peru and Brazil operated routinely during the reporting 
period; predictions were transmitted for approximately 140, 000 points on 3000 arcs, 
and nearly 26, 000 points of data wore received. 

Under the auspices of the NASA Kohoutek project, the Communications Center 
began total coverage in December in support of Comet Kohoutek. This round-the- 
clock operation will continue until February 1974. 

2.G Data Services 


The Data Services Division provided pointing predictions for four retro reflector 
satellites, BE-C, Geos 1, Geos 2, and Peole (this last satellite was covered during 
July and August only), for the laser sites at Mt. Hopldns, Brazil, Peru, South Africa, 
the NTU at Dionysos, and TAG. Gx'bital elements were routinely provided to Air Force 
Cambridge Research Laboratories, ONES, and the Institufe fiir Angewandte Geodasie, in 
Germany, for use in generating laser predictions. Camera predictions were generated 
and sent to 10 Baker-Nunn stations. Table 5 lists the satellites observed. 
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Tablo 5. Satellites tracked from 1 July through 31 December 1973 


salcimo 

Name 

Tracked on rcciuest from NASA 

IDOa ri3A 

lixiJorer 10 

i-oor, koa 

Geos 1* 

lOdfl 2A 

Geos 2* 

1070 loan 

Co.smos 362 rocket l)ody 

Tracked 

for geodesy and earth physics 

IflGl n51 

Midas d 

lOO'l C'lA 

nn-13 

lOOa 32A 

BK-C* 

IDflG 80A 

Geos 1 

liHlrt 2 A 

Geos 2 

1.070 10 OA 

Poole* 


Special requests 

1071 r.‘IA 

'I’hor burner 2 rocket 

1072 37A 

Molnlya 2 

1072 (59A 

Triad 

Tracked 

for long-period peitiirbations 

1050 al 

Vanguard 2 

1000 l2 

Echo 1 rocket body 

lOG-l 6.1 A 

BE-B 

1005 BOA 

Geos 1 

Trackcrl for atmospheric investigations 

1950 al 

Vanguard 2 

1900 el 

Explorer 8 

1063 53A 

Explorer 19 

I960 44A 

Explorer 32 

1008 06A 

Explorer 39 


Tracked by both lasers and cameras during 
this period. 
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Orbits and predictions were provided to the Baker-Nunn sites for satellite 1371 !34A 
In support of Dr. King-Hclo’s research and for satellite 1972 G9A (Triad) to help eval- 
uate the surface-force compensation system DISCOS as a future technique for gravity- 
field measurements. 

Local predictions were computed for SIcylab for public information. 

Final validation of SAO laser data was completed for Fiscal Year 1973 mul made 
available to the Analytic Satellite Geophysics Department for analysis. 

The film-control section received ajid cataloged 8476 observations on 7812 films 
from the Baker-Nunn stations. During this period, the practice of spooling three 
observations on one roll was started in order to conserve space and film cans. We 
received 840 films from the Air Force sites. 

Table 6 lists the 1520 precise reductions of satellite positions completed, wliich 
brings the number of all reductions as of 31 December 1973 to 242, 085. Of these, 

92 were provided by the group at the Instituto di Geodesia, Bologna, Italy, in coopera- 
tion with SAO. 

Photographs of 280 laser pulses were precisely reduced. 


Table 6. Reductions completed 1 July through 31 December 1973. 


Object 

Period 

Number 
of images 

Investigator 

Vanguard 2 

20 July to 20 August 1972 

52 

Jacchia 

Explorer 19 

20 July to 20 August 1972 

212 

Jacchia 

Explorer 32 

20 July to 20 August 1972 

103 

Jacchia 

Explorer 39 

20 July to 20 August 1972 

42 

Jacchia 

Explorer 19 

June through October 1972 

505 

Jacchia 

Explorer 39 

11-31 July 1968 

92 

Jacchia 

Pageos 

(selected 

observations) 

March .1967 to January 1970 

514 

Dr. I. Mueller 

Total 


1520 
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3. SATELLITE GEODESY AND GEOPHYSICS PROGRAMS 


3. 1 Introduction 


Durinp the period 1 July to 31, December 1973, wo began working on the 10 RTOPs 
assigned to SAO in support of EOPAP. These include studies of satellite dynamics to 
improve satellite ephemerides; use of existing satellite-tracldng data for determination 
of geophysical parameters, such as the gravity field, earth tides, polar motion, and 
station location; evaluation of SE HI; reduction and analysis of geodimeter measure- 
ments to determine local movements of the earth's crust; and studies to establish 
requirements for satellite-ti'acking equipment, network distribution, and satellite 
design. 

This period was one of transition from work supporting the NGSP to that for EOPAP. 
A number of NGSP activities were completed. The SAO contribution to the NGSP, 
based on SE m (Gaposclildn, 1973), was submitted to the American Geophysical Union 
and will be published as Chapter 9 in the NGSP Final Report. A comparative study of 
all NGSP contributions is nov/ under way. Preliminary results of this comparison 
indicate that geocentric coordinates have quoted accuracies optimistic by factors of 2 
to 5. Even so, independent estimates of coordinates agree with each other to within 
10 m, the NGSP goal. 

The scientific work related to individual RTOPs is discussed in the following 
subsections. 

3,2 Geophysical Data Base 

Under RTOP 3S9-01-04, Dr. M. R. Williamson and Gaposchkin have begun the 
develo,oment of computer software to manage a geophysical data base, which will be 
computer-accessible, through an interactive terminal if required. The data base is 
planned to have 10^ characters initially and to be expanded to 10^ or 10^^ characters. 
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The software being developed in FORTRAN for this data base Is a sot of user sub- 
routines; these routines will manage a data sot within the data base and are the same 
routines the agent who maintains the data base will use. Therefore, tlio user can 
create Ids own private data sets. The format and ordering of each data set is arbitrary, 
subject to the limit of about 100, 000 units of data in each set. Each unit can be 
arbitrarily large. 

The principal storage medium is a magnetic tape. A separate directory of infor- 
mation is maintained on disk. With this facility, it is easy to update and expand a data 
set or the data base and retain old data if desired. In fact, a user can determine 
quicldy what data are in the data base and in each data set without actually compiling 
subsets. Therefore, a user can retrieve a wide variety of data s.Jts, or subsets, 
depending on the contents of the data base. For example, for a data set of topographic 
heights, a subset can be requested containing all gravity and heat-flow measurements, 
where both are available, from places with a topographic height of > -1km. Finally, 
the facility was planned for the convenience of scientists; therefore, a worldng knowledge 
of FORTRAN programing is sufficient to use these subroutines. 

The first version of these subroutines is virtually complete. The subroutines are 
written in FORTRAN except for the random-access facility, and they can be transferred 
to another computer if a suitable random-access facility is available. A set of weather 
data are currently being incorporated into the data base. 

3. 3 Orbit- Computation Techniques 

RTOP 161-05-05, on orbit-computation techniques, encompasses two main activi- 
ties: 1) theory and computation of satellite motion and 2) development of minicom- 
puter software for field stations. 

3.3.1 Nongravitational perturbations 

The effects of direct solar radiation pressure and albedo radiation pressure on 
satellite ephemerides have been studied. While semianalytical developments for both 
these forces have been used for some years, the short-period effects have virtually 
been ignored uitiii now, owing to their small amplitude. Because the EOPAP program 
requires centimeter accuracies, we have been reexamining the adequacy of existing 
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theories. The analysis to date has been for spherical satellites having small or modest 
eccentricity, with particular attention given to Lagoos. The PageC'S satellite is being- 
used to test tWs formulation, as it has an area-to-mass ratio 10*^ times larger than 
Lageos's. 

Dr. D. A. Lautman has developed a theory for albedo perturbations that includes 
siioi-t-period effects. These formulas are now being incorporated into the orbit- 
computation program. Table 7 shows tlie approximate amplitudes of the albedo per- 
turbations of Pageos and Lageos expressed in linear units. The two orbits have the 
same sun-tiodc angle and the same argument of perigee. The largo magnitude of the 
along-treck secular change arises from the perturbation in the mean motion. A theory 
Is being developed to account for latitudinal variations in the albedo. 


Table 7. Theoretical calculation of albedo radiation-pressure perturbations. 


Parameter 

Pageos 

Lageos 

a 

1. 66 a 

1.9 

e 

0.06 

0.02 

I 

86f9 

90* 

a/m (cgs) 

136.0 

0.004 

Albedo 

0.4 

0.4 

Cross-track short-period perturbation 

10 m 

0. 05 cm 

Along-track short-period perturbation 

150 m 

1. 5 cm 

Cross-track secular perturbation 

100 m day"^ 

0. 25 cm day"^ 

Along-track secular perturbation 

3600 m day”^ 

25 cm day"^ 


The treatment of the direct effect of solar radiation has been revised by Dr. K. 
Aksnes to include an analytic formulation for both long- and short-period perturbations. 
This revision agrees well both with numerical integration and with the semianalytical 
formulation currently in use and is more efficient i,i terms of computer resources. 
These formulas will be tested in the orbit- computation program, and extension of the 
development is being considered for satellites with variable aspects. 
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3 , 3.2 Gravitational pertvjrbatlons 


Aksnes has completed a review of the analytical theory for gravitational pertux’ba 
tlons. A development of the earth's potential in spherical harmonics assumes the 
following fox’m: 



r 




00 I 

Pjj (sin 0) H- EE 

i=2 m=l 



P^^Jsin 0) cos(k-Kjj^^^ 


( 1 ) 

where the dominating zonal coefficients are 

Jg= 1.08X10"^ , Jg = -2.4XiO“® , J^= -1.7X10"® . (2) 

The higher order zonal coefficients J, and the tesseral coefficients J, are all of 
the second order in J or smaller. However, the perturbing effect of the correspond- 
ing harmonics on high satellites of geodetic interest diminishes rapidly with increasing 
values of I because of the quotient (R/r) , For high satellites, therefore, it is reason- 
able to assume that an orbit theory is adequate if it incorporates the first-order per- 
turbations due to all the known harmonics and the second-order perturbations arising 
from Jg, Jg, and J^. Furthermore, those tesseral harmonics that give rise to 
resonances in the motion of a given satellite also require a special nonlinear treatment. 

As a first step in testing the above supposition, Aksnes compared his second-order 
analjrtic theory (Jg, Jg, and only) with a theory of special perturbations (numerical 
integration). Complete to 0 (j 2 ), the analytic theory must account for all the following 
terms, most of which arise through interaction with JgS 

First order: Jg, Jg/Jg, J4/J2 

Second order; Jg, Jg, Jg/Jgj J4) J4A2* 

i 

The results of such a comparison are shown in Table 8 for Pageos (in orbit) and 
Lageos (to be launched). From column 5 , we see that for both satellites the maximum 
residual (= the error of the analytic theory) after 1 day of motion is about 20 cm, while 
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after 1 revolution {column G), It varies from 3 to 13 cm. During both intervals, the 
residuals vai’y roughly as sin 2u {u «= argument of latitude) and do not appear to increase 
rapidly with time. To find out which of the three harmonics contributed most to the 
error, two additional comparisons (last two columns) wex’o made for Lageos. Appar- 
ently, Jjj holds this distinction. 

Table 8. Comparison of theory and numerical Integration. 


SalollllD 

n (Mm) 

0 

1 

t 

A (cm dny"^) 

* 

A (cm rev" ') 

A (cm rev ) 

-1 * 

A (cm rev ) 

Pngeos 

10, !1 

0.0(5 

805 0 

10 

•1. 1 

- 

- 

Uigetifl 

10. 13 

0.02 

(50 

20 

13 

1.5 

0. 13 

Lpgoos 

10. 13 

0.02 

00 

21 

3.3 

1.2 

0.72 

Nolo A = 

"intoKriiled" 

mlnuB 

"tmalyllo" 

position. 





Jqi Jy) J.j Included, 
iind J I included, 
inuluded. 


The above is in agreement with the so-called d'Alembert characteristics, accord- 
ing to which the peiturbations in, for example, r must take the form 


6r = ^ e, i) e I ! sin I ^ I i ® Jjg (ju + kf) . (3) 

j,k 

When e is small, as in equation (3), terms with k = 0 must dominate. We conclude that 
an accuracy better than 20 cm requires the inclusion of third-order terms due to Jg in 
the analytic theory. For orbits with low eccentricity, only a few such terms are needed. 
An alternative to a third-order analytic development is to introduce some empirical 
terms of the form just shown. 
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Of course, it remains to be demonstrated what accuracies are achievable if the 
higher order zonal and the tessoral harmonics are also incorporated Into the theory 
(as well as into the numerical integration) in the manner already indicated. Most of 
these perturbations have already been included in our precision orbit-computation 
program. 

Satellite resonance has been studied by Ms. B, A. Romanowicz. Her develop- 
ment attempts to find a unified treatment for resonance that includes tlie resonant 
gravity-field coefficient and Jg by the Lie-Hori meihod. A first-order solution in 
elliptle integrals has been found. The formulation is being put into final form, and 
extension to second order is being investigated. 

3.3.3 Tides 

Dr. y. Kozai's revised theory for perturbations due to the sun and the moon 
(Kozai, 1973) has been incorporated into the precision orbit-computation program. 

Tliis theory includes parallactic and short-period terms and makes no assumptions 
about the motion of the sun or moon. The body tide and the associated phase angle 
are treated separately for the sun and the moon. At the same time, minor corrections 
to the perturbations due to the moving equator of the earth have been made. 

3.3.4 Numerical integration 

Williamson has converted the numerical-integration program previously used to 
verify the analytical formulas to an Adams method, a method far more efficient than 
that of Runga-Kutta-Gill. Investigation of the accuracy of this version continues. 

The program has been expanded to include an arbitrary gravity field expressed in 
spherical liarmonics. Modifications are in progress to include a surface-layer 
density model and the albedo radiation pressure in the force field. The integration 
program is being modified to be incorporated into the precision orbit-computation 
program in place of the analytical theory. This will provide a further comparative 
test with real data of the analytical treatment. 
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3. 3. 5 Mlnicomputor development 

Mr. J. R. Cherniack has continued directing the development of mlnicomputor 
software for tlie Data Genex’al Nova 1200 in support of field operations. Using orbital 
elements, the Nova computer field software has generated predictions that agree with 
tlie main computing-center predictions to 0?001. The programing required to aclUevc 
this agreement included the following; 

A. Assimilation of floating-point processors into FORTRAN. 

B. Incorporation of lunisolar perturbations. 

C. Incorporation of tesseral and zonal-harmonics perturbations by using a 
medium-sized field (16 X 16 plus resonant terms). 

D. Inclusion of estimated signal-strength computation. 

The programing was unusually difficult because of the limited memory size of the 
minicomputer, the lack of peripheral devices, and the absence of a decent operating 
system. These difficulties persuaded the research group not to use FORTRAN for 
subsequent field software. 

Additional software, which has been defined and is now being designed, will permit 
the following: 

A. Control of the laser mount in real time. 

B. Rejection of spurious data. 

C. Handling of pass summaries, report writing, etc. 

We decided to use FORTH (not to be confused with FORTRAN H or any other 
FORTRAN) for all the new software described above, A version of FORTH received 
from the National Radio Astronomy Observatory was upgraded in the following areas; 

A. Floating-point operations. 

B. Diagnostics and debugging, 

C. Addition of arrays. 

D. Mathematical library. 

E. Soft restart following power failures. 


19 



specifications and requirements for field minlcomputei’S have been established 
and documented {see SAO Proposal P 425-S-73, Supplement No. 1), Installations in 
the field will begin when authorized by NASA/OTDA. 

3. 4 S e lection of Fundamental Sites 

The analysis of ground-systems requirements and plans for a fundamental EOPAP 
nctwoi’k under RTOP 161-01-03 began. Assuming the number of stations and their 
distribution among the major plates are Imown, the criteria to be used for site selec- 
tion have been established and a particular site has been evaluated. Those criteria 
are discussed in the following sections, 

3.4. 1 Geological and geophysical environment 

The geological and geophysical environment is critical. A fundamental satellite- 
tracking site should be in a tectonically inactive region; that is, places where local 
deformation is occurring should be avoided. Local deformation can be studied more 
easily with other methods and will corrupt analysis of global motions. The principal 
types of information considered in choosing a site are as follows: 

A. Geological record. 

B. Seismicity of the area. 

C. Distance from the oceans. 

The geological record indicates the long-term (millions of years) evolution of a 
region, for example, by the existence of fault zones or uplift. Furthermore, old 
(Precambrian basement) rocks have recently been tectonically quiet. The seismicity 
of a region indicates the amount of strain release currently going on and, by implica- 
tion, the amount of deformation. However, anomalous earthquakes can occur in rela- 
tively quiet areas. Therefore, there is no foolproof method of selecting quiet sites, 
but the chances are best in an old craton or shield with no observed seismicity. 

Ocean tidal loading is a significant signal, which can be modeled only approxi- 
mately today. To minimize this effect, it is preferable for a site to be located several 
hundred kilometers from an ocean. 
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3. 4. 2 Atmospheric environment 

The principal considoi'ntion in selecting sites is the amount of cloud cover, as 
clouds arc the main reason for loss of data. Global cloud-cover maps are available 
(see, e, g. , Kravtsova, 1972; Landsberg et al . , J,9GG; Army-Air Force, 1942), and 
these can be used to identify good regions. More detailed information is then requested 
from local weather bureaus. In the future, the amount of relative humidity will dictate 
the accuracy of the refraction corrections, 

3.4.3 Orbital coverage 

Accurate orbit determination requh’es complete orbital coverage. Continuous 
observation is not necessary, but it is desirable that each section of the orbit be 
observed at least once in the interval of orbit computation. Orbital coverage depends 
on the distribution of stations, the satellite heiglti^ and the orbital inclination. 

Our station complement currently Is such that satellite orbits are covered ade- 
quately in the northern and equatorial latitudes and inadequately in the southern 
hemisphere. The critical parameter is the geocentric angle in the orbital plane 
delimiting the region of geometrical unobservability. The smaller the unobservable 
orbit-plane angle is, the better. This angle is a function of the southernmost station's 
latitude; we give its values in Table 9 for several satellites, for several choices of 
southernmost station, and for two assumptions about the minimum observable instru- 
ment elevation angle. 

At 5690 km, Lageos, with an inclination of 70°, will require a station latitude of 
25^7 for 15° elevation angle and of 29.°8 for 20° elevation angle. Therefore, complete 
Lageos orbital coverage is possible with any of the stations in Table 9, 

3.4.4 Communications and logistics 

Communications and logistics considerations require that stations be in inhabited 
regions and have communications. Further, it is desirable that the regions have com- 
mercial power and telephones and political stability. 
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Table 9. 

Southern unobservable orbit-plane angle in degrees. 
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3. 4. 5 Summary and recommendations 

Gaposchkin, Pearlman, and Dr. P. A, Mohr are collecting weather data and 
geological, tectonic, and seismic maps by regions. With more detailed data, they 
will investigate the number and distribution of the tectonic plates among a complex of 
stations. 

The necessity of moving the laser station 7902 out of South Africa provides an 
opportunity to apply these considerations. Because of its advantageous latitude, 
Australia has been studied in some detail, and two sites are documented here. Geo- 
logical and seismic maps of Australia are presented in Figures 1 and 2 (McElhinny, 
1973). The first site studied is in Orroral Valley, where there is an existing NASA 
facility and a lunar laser system. However, Orroral Valley is not suitable as a long- 
term EOPAP station, because it is directly on a seismic zone in the midst of a Paleo- 
zoic erogenic region. Since South Africa is considered to be a good station for observ- 
ing, we compare its cloud-cover data (Figure 3) and its data loss due to cloud cover 
(Figure 4) with those for Orroral Valley (Figure 5) (Morgan and Miller, 1973). It is 
clear that there would be a significant loss of data due to cloud cover in Orroral Valley, 
yet from Table 9, we can see that the orbital coverage for that region is an improve- 
ment over that in South Africa. 
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Figure 1, Seismic map of Australia. 
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Figure 2. Geologic map of Australia, 
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Figuro 6, Mean cloudiness, Orroral Valley. 


The region around Kalgoorlie in southwest Australia is the second site studied. 
Figures 6 and 7 show the local cloud cover for a typical Eastern seacoast area of 
Australia near Orroral Valley (Sydney) and for the southwestern inland region near 
Kalgoorile (Lawlers). The weather is far superior in both these areas to that in South 
Africa and in Orroral Valley. From geological and seismic maps, it is apparent that 
tliis region is located on one of the oldest Precambrlan basement rocks Imown («3 billion 
years) and has no seismicity. Because of gold mining in that region, good communica- 
tions exist there. In addition, the mines provide the opportunity to establish stations 
to monitor earth tides. Finally, the orbital coverage offered by this site is acceptable. 
Gaposclilvin, Pearlman, and Mohr are obtaining detailed geological and tectonic maps, 
other geophysical information, and local weather data of this region. Tentatively, 
however, it appears that locating a laser tracking station in the Kalgoorlie region would 
be ideal for EOPAP from the scientific point of view. 
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3.5 Gravity Field 


Work on the earth's gravity field, under RTOP 369-02-03, has been concerned 
with the analysis of long- term effects, Gaposohkin and Ms. G. Mendes have recom- 
puted a long series of orbits for Geos 1 and Geos 2 by using Baker- Nunn data. These 
series have been examined for the following; 

A. To test the lunar, solar, tidal-radiation-pressure, reference-system, and 
long-period gravitational -perturbation theories. 

B. To evaluate the internal coherence of mean elements, 

C. To pi'ovide reduced mean elements for the determination of those harmonics 
in the gravity field that give rise to long-period perturbations. 

The internal eonsistency of mean elements depends on the averaging time - the 

number of days covered by each orbit. With Baker-Nunn data, the low density of 

“5 

observations precludes calculating orbits coherent to 10 with 2 days of data. A 

-5 -6 

4-day span of data will give 10 to 10 , and except for unusually poor distribution, 
an 8 -day span will give a, e, and I consistent to lO"*^, 

The reference-system correction, tides, and parallactic terms in the lunar and 
solar perturbations all have the same amplitude and frequency. By using Love number 
kg compensated for the ocean tide, these perturbations seem to be treated properly. 

Residual long-term variations in a, e, and I remain unexplained. The series of 
data is being extended to a full year. These residuals could be due to the as-yet- 
unmodeled albedo radiation (RTOP 161-05-05) or to the anisotropic reflectivity and 
variable cross section of the satellite. Both possibilities are under investigation. 

From tills analysis, an unexplained long-period effect in the eccentricity related 
to the odd zonal harmonics has been observed. This effect cannot be due to tides or 
reference systems, as these do not influence the eccentricity. The neglected second- 
order terms in gravitation theory are currently being examined 

(RTOP 161-05-05) as a possible explanation. 



3. 6 Station -Coordinate Determination 


Under RTOP 3G9-02-01, two activities concerning tectonic plate and fault motion 
have been started; 

A. The use of data from the BAGEX program to determine the geocentric coordi- 
nates of individual stations, 

B. The design and development of a general simulation program, 

ISAGEX laser data on Geos 1 and Geos 2 have beej'i used to determine the station 
coordinates of Dakar (7820) and Haute Provence (7809), The coordinates were deter- 
mined in the same computation with the orbit in the precision orbit-determination 
calculation. 

The ISAGEX data for Dakar did not give a satisfactory solution for the station 
coordinates in SE in. The coordinates determined in SE m disagreed with the height 
of mean sea level and with the coordinates from the nearby BC-4 station by more 
than 40 m. We thought this subsequent investigation would provide acceptable coordi- 
nates and indicate what problems existed with the data. The study was also extended 
to Haute Provence because we liave reasonable confidence in the coordinates and laser 
data from that station. The results using current data, geodetic parameters, and 
analytical models with nine arcs of data can be summarized as follows; 

A. Only single-pass determinations were possible from 7820 owing to a paucity 
of data. 

B, The mean of these independent determinations agreed with the "best" global 
solution to within 15 m for station 7809 and 23 m for 7820, 

C, The independent determinations of coordinates differed from the "best" global 
solution by as much as 121 m for station 7820 and 44 m for 7809. 

D. The amplification of errors is evident from the study. The overall accuracy 
of thes'T 4-day ISAGEX orbits is about 7 m. Single-pass station determinations used 
for data screening and preliminary simulations consistently give a factor-of-10 ampli- 
fication (RTOP 161-05-02). 

It is evident that adding more arcs to the average will result in a better global 
solution. 
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The use of existing data has raised a number of questions that Gaposchldn and 
Mr. J. Latimer are now examining through a simulation study. For example, they 
plan to study the effects on station-coordinate determination of known gravity-field 
model errors, the distribution and errors of station coordinates, the number of sta- 
tions and Icnown eri’ors In their data (such as refraction, retrorefloctor-array distor- 
tion, and signal strength), polar motion, mean orbital elements, and the gravity field. 
The first use of the simulation will be to establish the number of stations needed to 
measure plate motions. The program being instructed will employ the precision 
geodesy program for all analyses. The simulation consists of defining the universe, 
gravity-field station coordinates, tides, drag, radiation pressure, and system errors 
and generating corresponding data. These data can be corrupted (e. g., by loss due to 
weather, system biases, poor station distribution) and used to compute geophysical 
and geodetic parameters. Shortly, it will be possible to perturb each system param- 
eter separately. The influence of critical errors - such as station distribution, orbit 
configuration, gravity-field uncertainty, and refraction errors - can )3C determined in 
this way. 

3. 7 Surface Surveys 

Under RTOP 161-05-06, Mohr and Mr. A. Girnius have made further analyses 
of all the Ethiopian geodimeter surveys to date (Mohr, 1974). The precision of the 
survey lines has been improved by the derivation and application of an empirical 
formula for ground- radiation effects, to be applied to the atmospheric-density correc- 
tion. This precision is now ±5 to 6 mm. However, under stable atmospheric condi- 
tions, the precision may improve to ±3 mm, and, thus, crustal deformation of this 
magnitude should be detectable. 

The new analysis reveals progressive crustal extension in the western sector of 
the northern network at a rate of 1.0 cm yr”^. This extension acts perpendicular to 
the rift trend, as expected from geological considerations. The central Ethiopian rift 
network has suffered progressive deformation at a rate of about 0. 5 cm yr \ again 
perpendicular to the rift trend. The southern network has not been significantly 
deformed. The northern network also reveals possible rift-trend zones of no extension 
and even of compression; the application of a homogeneous strain field, commonly 
assumed in analyses of geodimeter surveys of the San Andreas fault region, California, 
appears unjustifiable for a grid-faulted rift valley. 


29 



The deformation rates oV)ia5iwed from the Ethiopian geodimetor surveys are not 
inconsistent with estimates Mohr has made from measurements of the magnitude of 
rift faulting and from radiometric ages of the rocks cut by the faulting (Mohr, 1973). 
Mohr concludes that the rate of crustal extension at the African rift plate boundary 
has quickened during the Quaternary, and tliat the zone affected by extension has 
progressively narrowed to its present 5- to 10-km width within the rift floor. 

Worldng with a survey team from Dartmouth College and the University of Iceland, 
Mohr has remeasured geodlmeter lines across the Eastern active zone of Iceland. 
Previous surveys revealed an extension rate of 2 cm yr"^ for this zone, in precise 
agreement with plate -tectonic studies of the North Atlantic (Decker et al ., 1971). 

But the 1973 resurvey showed a more complex picture: Sectors that had extended 
from 1967 to 1970 tended to shorten again between 1970 and 1973. This indicates that 
rift deformation on a scale of years only is unlikely to match integrated plate-tectonic 
rates. It is lioped that in this new field, precise geodetic surveys will suggest the 
manner and cause of plate separation. 

3. 8 Polar Motion and Earth Rotation 


The study of polar motion and earth rotation, RTOP 161-02-02, is concerned with 
determining the accuracy of mean elements (RTOPs 369-02-03 and 161-05-05) and the 
preparation of EPSOC data for analysis (see Section 2. 1). 

By using 8-day averages of Baker-Nunn camera data, the mean inclination of 
Geos 1 and Geos 2 can be computed with a consistency of 10 . The EPSOC data for 

Fiscal Ye i. 1973 are now preprocessed and can be used to determine mean elements. 
The success in improving the accuracy of mean elements in general and the mean 
inclination in particular depends on the proper treatment of perturbing effects and on 
the accuracy of tracking data. Of particular importance now are the effects of photon 
pressure and body and ocean tides (RTOPs 369-02-03 and 161-05-05). The EPSOC 
data will be more consistent and of greater accuracy and reliability than any laser data 
used before. Preliminary orbits are now being computed. 

The theoretical model of polar motion by Dr. G. Colombo is progressing. This 
deterministic model is based on nonlinear coupling between the core and mantle to 
provide a means for transmitting the excitation energy for the Chandler motion. The 
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effects that earthquakes have on tlie Chandler motion are currently being investigated, 
not as a source of excitation energy, but as a source of random changes in the inertia 
tensor. The speetral characteristic of the Chandler motion excited by core-mantle 
coupling and perturbed by earthquakes is being compared with the observed spectrum 
obtained by the maximum -entropy method. The preliminary results are very 
encouraging. 

3. 9 Laser Techniques 

Under RTOP 161-05-02, the first draft of a detailed description of the techniques 
used by SAG to compute retroreflector-array transfer functions has been completed 
by Arnold and will sot>72 be published as an SAG Special Report, The text includes 
methods to account for cube-comer size, shape, and optical properties, as well as 
discussions on the influences of Interference and diffraction, A copy of the draft was 
furnished to GSFC for their review and use. Preliminary retroreflector-array transfer 
functions were calculated for Geos C and Starlette during Fiscal Year 1973. The final 
ti’ansfer functions will be computed as soon as complete spacecraft designs have been 
provided by the Applied Physics Laboratory and ONES. 

Several analyses have been carried out by Arnold and Pearlman in support of the 
Lageos program. Simulations based on the retroreflector-array transfer functions 
were used to investigate a number of spacecraft designs and orbital parameters to 
ensure that Lageos will satisfy program requirements and yet be compatible with 
projected laser capabilities. Simulations and field experiments were also used to 
study the ability of the Baker-Nunn camera to track Lageos with a number of different 
spacecraft and orbital parameters. Tracking with the Baker-Nunn camera will be 
required for routine orbital maintenance and for laser pointing predictions. The results 
of the SAG analyses have been furnished to Dr. J. Siry and Mr. R. Spenser of NASA 
Headquarters and to Mr. J. Randell of NASA/MSFC, 

System simulations were performed on models of the SAG laser systems to eval- 
uate specific designs for laser upgrading for Geos C. These simulations indicated 
that pulse digitizing techniques would reduce bias errors introduced by the present 
leading edge, fixed-threshold detection techniques to less than 10 cm. These analyses 
were basic support in the SAG proposal to GTDA for laser upgrading. 



In accordance with a plan to achieve increased range accuracy for Geos C meas- 
urements, a study of the systematic and random errors of the SAG laser systems in 
their current configuration has been made by Mr. C. G, Lehr and Pearlman. In their 
present configuration, the SAG systems determine range from readings of a time- 
interval counter that is started and stopped at threshold points on the leading edges of 
the transmitted and received pulses, respectively. Since the counter's resolution 
(0. i or 1 nsec) is considerably less than the laser's pulse width (20 to 25 nsec), the 
counter readings are sensitive to the variations in the size and shape of the returned 
pulses. A random error of 75 cm is typical; but of greater concern is the systematic 
error, which is insensitive to statistical averaging. Figure 8 shows how this systematic 
error varied during a transit of Geos 1 at Arequipa, Peru, The correction for the 
systematic error is largest near culmination of the satellite, because at this time the 
signal's amplitude is greatest and consequently the counter's threshold level is rela- 
tively farther down the pulse's leading edge. The correction curve of Figure 8 was 
obtained by measuring oscilloscope photographs of return pulses. Since these photo- 
graphs capture the amplitude variations within the pulse, they can be used to compute 
the time between threshold crossing and pulse center. Both the effects of intensity 
variations from pulse to pulse and the quantal irregularities within pulses are accounted 
for by this procedure. 
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Figure 8. Example of systematic range correction, taken on Geos 1 transit^ 9 June 
1973 from Arequipa, Pem. 
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The photographic method is not efficient at a high data rate but is very useful for 
system analysis on selected transits. Investigations based on pulse photography were 
used in the development of SAO’s plan for laser system upgrading for Geos C. A 
waveform digitizer, soon to be Installed at the Mt. Iloplcins observing station, will 
take data automatically and correct for irregularities in pulse shape and size. 

In a program to understand the opei’atlon of individual system components and to 
identify any unsuspected error, "s, measured and calculated values of the numbers of 
photoelectrons were compared for return pulses obtained from satellites and fixed 
reflectors. It was found that the photomultiplier tubes show saturation effects at 
about 1000 photoelectrons. As part of this investigation, receiver efficiencies were 
measured by using stellar photometry. The values (expressed as photo electrons 
emitted by the photomultiplier tube divided by photons entering the receiver aperture) 
are listed in Table 10. 


Table 10. Receiver efficiency measurements. 


Station number 

Location 

Receiver efficiency (%) 

7902 

South Africa 

1.9 

7907 

Peru 

0.6 

7921 

Arizona 

1.1 

7929 

Brazil 

3.2 


3. 10 Other RTOPs 

Work on analytical models of earth motions, RTOP 161-02-07, will begin in the 
next few months. Satellite self-tracking by use of altimetry data, RTOP 161-03-07, 
will be combined with the Geos C sea-surface topography program. Work on the 
Geos C program is expected to begin during the last part of Fiscal Year 1974, by which 
time the NASA/Wallops contract for Geos C will have been received by SAO. 
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4. A TMOSPHERIC RESEARCH 


Jacchla’s attempts to model the variations in thermospheric composition as 
observed by mass spectrometers on the OGO 6 and ESRO 4 satellites and to reconcile 
them with the dopplor temperature data from OGO 6 and total densities from satellite 
drag have scored remarkable success (Jacchia, 1973). Two simple formulas have 
been found to represent the observed seasonal-latitudinal variations in composition and 
the diunial variation of the individual atmospheric constituents at any given location 
and height. Each of the equations contains a single numerical constant. 

Wox’k Is in progress to determine with greater accuracy the numerical value of the 
constants. Also, attempts are being made to find a way to represent the global varia- 
tions of temperature and composition with geomagnetic activity. In the latter project, 
a collaboration is envisaged with a group at the Physics Department of Bonn University, 
West Germany, which is analyzing the ESRO 4 composition data. 

A study of the effects of radiation pressure and atmospheric drag on the balloon 
satellite 1963 30D has been completed (Slowey, 1974). Slowey has found that by using 
a single value of the reflection parameter, it is possible to model the perturbations due 
to direct radiation pressure in all the elements. This is contrary to the earlier findings 
of other workers and is thought to be due primarily to the more accurate treatment 
given to terrestrial radiation pressure. The reflection parameter was found to change 
very little from the beginning to the end of the 7-yr lifetime. A very small increase 
was detected, however, which may be related to a slight change in the character of the 
surface. Another interesting result was the success achieved in modeling the "anom- 
alous" perturbations that now seem clearly associated with the effect of direct radiation 
pressure on the aspherical figure of the satellite. These perturbations are quite 
noticeable in the orbital acceleration and were also clearly detected in the eccentricity. 

A rotating component of force, transverse to the solar direction, was found to reproduce 
the effects extremely well. The magnitude of this force was only about 3% of the total 
force due to direct radiation pressure. It, too, remained essentially constant over 
the lifetime, although the period of rotation changed considerably. 
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Terrestrial radiation pressure can make a significant contribution to the total 
acceleration of satellites at greater heights in the earth’s atmosphere. The success 
attained using the Tiros model of the earth's albedo in the case of 1903 SOD led the 
atmospheric research group to apply the model to the orbits of other liigh satellites 
the determination of atmospheric drag. At the same time, a more precise assess- 
ment of the effects of direct radiation pressure on these satellites was made. As a 
result, a catalog giving more accurate densities from five balloon satellites will bo 
published in tlxe near future. Work has now begun to utilize those results in new 
studies of the hydrogen and helium concentrations and of the oxygen interface. Much 
of tills, of course, is related to the determination of constants in Jaccliia's new model 
for the seasonal -latitudinal variations of all atmospheric constitu >nts, a project that 
is just beginning. 

Slowey has begun work on tlie determination of atmospheric rotation from the 
secular decrease in the inclination of satellite orbits. Most of the necessary program- 
ing has been completed, and the first results were being obtained at the end of the 
reporting period. The initial work uses the periods of orbital decay for four balloon 
satellites. Results indicate that significant values of the rotation can be achieved 
from 2 to 3 months of data in these cases. The values obtained apply to heights where 
few other data are available, and these values are quite important to an understanding 
of the rotation. If the method is successful, it is planned to use six or seven ordinary 
satellites in lower orbits to obtain mean values of the rotation. 

Normally, work would have continued on the determination of atmospheric densities 
from the drag on six or seven artificial satellites with perigee heights ranging from 
130 to 800 km. However, funding restraints, particularly in the areas of personnel and 
computing time, greatly impeded this study during the reporting period. It is hoped 
that we can find a way to sustain the density deteianlnations on a reasonably up-to-date 
basis. Assuming that this can be accomplished, the upper atmospheric research group 
would like to perform correlative studies by using extreme -ultraviolet data from the 
Atmospheric Explorer satellite and global density data from drag on a variety of satel- 
lites at different heights and under different conditions. 
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